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ABSTRACT 



Type la supernovae (SNe la), produced by the thermonuclear explosion of 
white dwarf (WD) stars, are used here to derive extragalactic distances and an 
estimate of the Hubble constant from their emission signatures at late phases 
{Nebular SNe la Method^ NSM). The method, first developed in Ruiz-Lapuente 
& Lucy (1992), makes use here of an improved modeling of the forbidden line 
emission at late phases. Hydrodynamic models of the explosion of WDs of 
different masses, both sub-Chandrasekhar and Chandrasekhar, provide the 
basis for comparison with the observations. It is shown that it is possible to 
probe the overall density structure of the ejecta and the mass of the exploding 
WD by the effect that the electron density profile has in shaping the forbidden 
line emission of the iron ions, and that a robust diagnostic of the mass of the 
exploding WD can be obtained. Cosmic distance scale can thus be related to 
basic diagnostics of excitation of iron lines. Once the most adequate model 
is selected, comparison of the predicted line emission at these phases with 
the observed spectra gives an internal estimate of both the reddening and the 
distance to the SNe la. The results presented here favor denser models than 
those corresponding to sub-Chandrasekhar explosions. From a sample of seven 
SNe la in Leo, Virgo, Fornax and beyond, a value of the Hubble constant 
Hq = 68 lb 6 (stat) lb 7 (syst) km Mpc~^ is derived. The depth of the Virgo 
cluster is found to be large, ranging from 13 to 23 Mpc at least. If NGC 4526 
traces well the core of the Virgo Cluster, then the latter is located at 16 ±2 Mpc. 
The galaxy NGC 3267 in Leo appears to be located at 9.8 ±1.5 Mpc. 



Subject headings: cosmology: cosmic distance scale — supernovae: general 
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1. Introduction 



The Hubble constant is a key to the model of our Universe. The long-standing debate 
on its value has recently been rekindled by the use of the HST to observe extragalactic 
Cepheids (Freedman et al. 1994; Sandage et al. 1994; Tanvir et al. 1995). When comparing 
results from different methods to determine extragalactic distances, the use of the luminosity 
of Type la supernovae (SNe la) at maximum as a "standard candle" has been at odds with 
other widely used distance indicators (Branch & Miller 1993; Tammann & Sandage 1995). 
Only recently it has been clearly confirmed the correlation of brightness at maximum with 
rate of decline of the light curve of SNe la, earlier suggested by Pskovskii (1984). The use 
of maximum brightness of SNe la has to be reconsidered on the light of this correlation 
(Phillips 1993; Hamuy et al. 1995; Riess, Press, & Kirshner 1995). All empirical uses of 
SNe la luminosities require a zero-point calibration via Cepheid distances to a galaxy host 
of a SNe la. A different approach, which gives the distances to SNe la in a single step, 
is to use theoretical light curves (Arnett, Branch, & Wheeler 1985; Hofiich & Khokhlov 
1995). However, the problem there is the existing limitations in discriminating among the 
underlying explosion models. 

Here we present also a "theoretical" approach (no need of any zero-point calibration 
and based on the assessment of models), but it differs from the light curve approach: to 
determine distances to SNe la we use the information contained in the emission of the 
SN ejecta at late phases. As in other astrophysical contexts, forbidden emission can be 
exploited here, since such emission of iron ions (Fe"'" and Fe"'""'") traces very well the electron 
density (ue) and electron temperature (Te) of the supernova "nebula". This fine sensor 
of density profile probes, as we will show, the mass of the exploding white dwarf. At 
these phases, the ejecta are kept warm by thermalization of the 7-rays and positrons (e"*") 
from the decay of ^^Co (coming from ^^Ni) into ^^Fe. This is the powerhouse which sets 
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the temperature of the ejecta. As in other "nebular" environments, the emissivity of the 
different forbidden transitions of Fe"*" keeps track of the electron density profile whereas 
both Fe"'" and Fe give information on Tg. The large number of forbidden lines over the 
wavelength range at which SNe la are observed gives, in addition, enough information to 
obtain an internal estimate of the reddening. The basis of the method has been presented 
in Ruiz-Lapuente & Lucy (f992. Paper I), and in later work the treatment of radiation 
processes has been refined (see Ruiz-Lapuente f995, R95, for a review). The input for 
the present study is provided by different models of the explosions, some of which involve 
multidimensional hydrodynamic calculations. 

The explosion mechanism of SN la is currently being reanalyzed on basis to more 
consistent physical treatments than the earlier phenomenological approaches. Two classes 
of models for the thermonuclear explosion of WDs as SNe la are now competing: (f ) central 
ignition of a C+0 WD close to the Chandrasekhar mass, and (2) detonation of a C+0 
WD below the Chandrasekhar mass (sub-Chandrasekhar WD) induced by detonation of a 
helium layer accumulated on the surface. The issue is relevant to the determination of i^o, 
since as addressed in Livne & Arnett (f995), if sub-Chandrasekhar WDs do explode, the 
distance scale derived from SNe la should be revised: light curves would place SNe la at 
shorter distances. The longer distance scale is linked to the explosion of C+0 WDs reaching 
the Chandrasekhar mass (Arnett, Branch & Wheeler f985). An important goal of current 
discussions on explosion models is thus to clarify which distance scale is actually favored by 
SNe la. In the present work, by addressing the distance scale problem through modeling of 
the late phases of SNe la, we equally obtain insight on the masses of the exploding WDs. 
In this Letter we apply the Nebular SNe la Method (NSM) — or Co-Powered Nebular 
Candle Method (more briefiy Nebular Candle Method^ NCM), as the reader might prefer — , 
to a sample of SNe la, well observed in their nebular phases. These SNe la are in galaxies 
of the Virgo Cluster, Fornax, and Leo Group. Our method has the advantage that it does 
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not depend on any zero-point calibration of the luminosities, and it provides individual 
distances to each SN la without assuming them to be "standard candles". It uses the "fine 
sensitivity" of nebular diagnostics for ejecta which, at the phases considered (about 8-10 
months after explosion) do not show any evidence of dumpiness or asymmetry in their 
spectral lines. A major advantage of our method is the possibility of directly probing the 
electron density profile ne(r) and the mass of the exploding WD. This possibility arises 
from the fact that the electron density at the center of the ejecta ne(0) is directly linked to 
the mass of the exploding WD, and the forbidden emission of the Fe ions is a tracer of the 
density profile of the ejecta down to the center. 



2. Physical model 



The basic idea of our approach, presented in Paper I, has later been developed with the 
use of realistic density profiles for the underlying explosions and by significantly improving 
the physical treatment of the radiation transport. In the present calculations we include 
a numerical treatment of the escape of positrons, which takes into account their energy 
spectrum and follows their path throughout the ejecta for different configurations of the 
magnetic field: a turbulent magnetic field and a radially combed out one. In this way we 
estimate the amount of positron energy deposited in the different layers of the ejecta. A 
detailed treatment of 7-ray transport is also done, but positron deposition of kinetic energy 
is dominant at the times considered. Nonthermal processes resulting from the various 
interactions of the 7-ray photons and positrons from ^^Co decay with the electrons and ions 
of the supernova plasma enter in the ionization balance. The ionization state is calculated 
with due account of photoionization arising in the recombination of the different ions. 
The rate equations are solved for forty most abundant ions. The electron temperature is 
self-consistently calculated from the balance between cooling by forbidden line emission and 



- 6 - 



heating by deposition of energy by 7-rays and positrons (Ruiz-Lapuente, Kirshner, Phillips, 
et al. 1995, RKP95). We include all fine structure transitions up to the term a^I for Fell 
and up to the term "^D for Fe III, for which accurate coUisional data from the Opacity 
Project exist (Pradhan & Berrington 1993; Bautista & Pradhan 1995). A full description 
of the physics can be found elsewhere (R95), and a complete discussion will be done in 
subsequent papers. The key diagnostics on density, however, are beyond the uncertainties 
on either atomic data or ion models. A basic "form-factor" discriminating between models 
is linked to excitation, where uncertainties of atomic data are well delimited. For the study 
of the density diagnostics, rate-equation decoupled from the rest of the treatment is used 
for complex Fe"*" and Fe models (54 fine-structure energy levels are included for Fe"*" and 
34 for Fe"*""*") to study the variation of emissivity with electron density. 

The outcome of central ignition of C+0 WDs with masses close to the Chandrasekhar 
mass has been calculated with one-dimensional hydrodynamic codes. The resulting models, 
such as the popular W7 (Nomoto, Thielemann, & Yokoi 1984), or DD4 (Woosley & 
Weaver 1995), although lacking a consistent treatment of thermonuclear fiame instabilities 
and propagation, have achieved reasonable success in their comparison with observations 
(Branch et al. 1985; Fastman & Pinto 1993; RKP95). In the alternative mechanism for SNe 
la, detonation of sub-Chandrasekhar C+O WDs from compression by helium detonation 
in an outer shell, WDs over a broad mass range would explode (masses from 0.6 to 1.2 
M0). This has been explored by different authors (Livne & Glasner 1991; Ruiz-Lapuente 
et al. 1993; Woosley & Weaver 1994; Nomoto 1995). Two-dimensional calculations coupled 
with nuclear reaction networks have recently been completed by Livne & Arnett (1995). 
The corresponding spectra are investigated in the present work (see also R95, and work in 
preparation with F. Livne). The set of models investigated here represent the whole class of 
so-far discussed mechanisms and progenitor masses for SNe la (WD masses ranging from 
0.6 to 1.38 M0). The observational sample consists of well-documented SNe la, for which 
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good, well-calibrated nebular spectra exist. Estimates of the errors in the flux scale are 
available (see, for critical observational issues Ruiz-Lapuente & Filippenko f994). 

3. WD mass and distance scale 

The density of the SNe la ejecta and Hq.- Diagnostics on total mass and electron density 
profile can provide a crucial test for the current debate on SNe la models. In late phases, 
the population of the energy levels of the ions is out of LTE across a significant part of the 
ejecta, because the density of the electrons responsible for coUisionally exciting the lines is 
lower than the critical density for the corresponding transitions. A clear diagnostic of low 
Ue comes from the drop of the emission at A5200 A, A4300 A, and A5000 A. These emissions 
are due to the Fe"'" a'*F-b'*P and a'*F-a'*H transitions, whose lower energy terms can be 
significantly depopulated if Ue is low. This makes the coUisional excitation rates decrease 
significantly. In this respect, an interesting ratio is that of a'*F-b'*P to a'*F-a'*P, with 
emissions at A5262 and A8617 respectively. As illustrated in Figure 1, this ratio becomes a 
diagnostic of mass. The infrared transitions provide a fine diagnostic on Ue too: a suitable 
one is the ratio of the [Fe II]A8617 to the [Fe II] 1.257 fim line. We restrict ourselves to 
the optical range, however, where most of the observations are being done. The emission 
ratios of different transitions of the same ion provide a robust test that reaches well beyond 
uncertainties in ionization balance, and give an approach to Ue which is independent from 
the other probes (ne(0) is basically a function of the mass: the number density at the 
center, n(0), varies by a factor larger than 10 between a 0.6 and a 1.38 WD whereas the 
electron fraction Xe is 1.5-2 along the luminous SN ejecta). In addition to the diagnostic 
on Ue, the emission of Fe"*" also bears interesting diagnostics of electron temperature (the 
ratio of the [Fe II] A7157, in the R-band, to those in the B-band, [Fe II] A4416, or in the 
I-band, such as the [Fe II] A8617 emission). 
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It is found in our analysis that, at the current stage of development of models, the 
sub-Chandrasekhar ones appear to fall outside the region of the (ue, Te) plane favored by 
the observations of "normal" SNe la. That produces discrepancies between the predicted 
relative luminosities of the lines and those displayed by the SN. Figure 2 (top panel) shows 
the calculated spectrum for W7 at 300 days after explosion, compared with the spectrum 
of SN 1994D, a normal SN la. A comparison of the same SN with the sub-Chandrasekhar 
model showing the closest resemblance to the observations (explosion of a 0.8 C+0 
WD after accretion of 0.17 of He) is done in the bottom panel. The picture where 
the emission of some Fe"'" transitions drops to zero, in SNe la from low-mass WDs, is at 
odds with the observations (Figure 2, bottom panel). Such a test allows us to conclude on 
the model density (massive or less massive WD) quite independently from the accuracy 
in the ionization balance. Our generic concern, valid for the whole range of He-induced 
detonation models, comes from the effects that the lower densities and, in many cases, the 
smaller ^^Ni mass, produce in the forbidden emission at late phases. The factor by which 
the distance scale would be shifted by adopting the best sub-Chandrasekhar models instead 
of the Chandrasekhar ones is 1.1. The Hubble constant would thus be centered at 75 km 
s~^ Mpc~^, instead of being centered at 68 km s~^ Mpc~^ as for the Chandrasekhar model. 
Such a shift, however, is disfavored by "normal" SNe la. In addition, the variety of models 
that He-induced detonations naturally predict does not seem to have a clear counterpart 
in the observations: explosions producing 0.3 of ^^Ni in a low-density WD, or those 
with I.IM0 of ^^Ni moving at high speeds have not yet been identified. Binary evolution 
effects might, however, play a role in selecting a particular, narrow range of WD masses for 
exploding by He-detonations (Canal, Ruiz-Lapuente, & Burkert 1996, CRB96). 

From the better agreement of "normal" SNe la with Chandrasekhar models, we obtain 
distances to our sample of SNe la, together with internal estimates of reddening. Table 
1 gives the distances and reddenings derived for the different SNe la, and the absolute 
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magnitude at maximum M^^^ deduced from these estimates. The value of the Hubble 
constant obtained by this method and from this sample is Hq = 68 it 6[stat) km Mpc~^ 
in concordance with the M^^^ implied for the SNe la. The error bar quoted here includes 
the effect of errors in absolute flux calibration and wavelength-dependent fluxes in 
the observations, and uncertainties in the internal estimate of reddening. In addition, 
we account for a 10% systematic error in Hq coming from the method (i.e. modeling 
uncertainties attached to photoionization and determination of the electron temperature). 
Our estimate of the Hubble constant results in Hq = 68 it 6[stat) it 7[syst) km Mpc~^ . 
( M^„^ ) for "normal" SNe la is -19.2 ± 0.3 (SN I99IT is about 0.3 mag brighter: M^„^= 
-19.53 ± 0.22). 

Extragalactic distances and SNe la.- The Nebular SNe la Method provides a new way to 
determine distances connected to our understanding of SNe la. The extragalactic distance 
scale given by this method, at its present stage, already indicates the considerable depth 
of the Virgo cluster, which should extend at least from fi = 30.6 to 31.8 (13 to 23 Mpc), 
since our subsample of SNe la located in Virgo covers such a range of distances. Fornax is, 
according to our result, at a distance equivalent to that of the W Cloud, in the background 
of the Virgo cluster. This distance is likely to be overestimated in Table I: SN I992A might 
be slightly less dense than what is predicted by the Chandrasekhar models considered here. 
Which is the distance to the Virgo core? Tammann & Sandage (1995) argue that NGC 4526 
being an early-type galaxy, it is likely to trace well the distance to the core of the Virgo 
cluster. Given our estimate of the distance to NGC 4526, the distance to the core would be 
of the order of 16 i: 2 Mpc. 

SNe la have been for a long time the stronghold of the long extragalactic distance 
scale. However, when using the expression "long" or "short" for astronomical distances, one 
should keep in mind the depth of the Virgo Cluster, and the uncertainties arising when a 
Hubble flow velocity has to be assigned to the above galaxies. This uncertainty is inherent 
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to any determination. Given this situation, it seems necessary to compare distances and 
absolute magnitudes for SNe la implied by different methods. 

How does the mass of the WD relate to the value of the Hubble constant? Values lower 
than 50 km s~^ Mpc~^ are excluded by the Chandrasekhar limit for the explosion of a C+0 
WD. The ejecta can not trap more energy, and the highest electron density attainable to 
populate the forbidden transitions is also limited by the Chandrasekhar mass: Ue (0) can 
not be larger than a few times 10^ cm~"^ at 300 days. The luminosity and the emissivity in 
forbidden lines is thus constrained by this requirement. From observations, it is also clear 
that the proper range of electron densities is below that limit: Fe"'" forbidden transitions 
would appear too strong for a density higher than 10^ cm at 300 days. Thus, Hq <^ 50 
km s~^ Mpc~^ is only possible if SNe la would come from objects more massive than the 
Chandrasekhar mass. An upper limit to Hq is equally set by the disagreement of explosions 
of light WDs with the bulk of SNe la observed (exception to this disagreement is only found 
in some very subluminous SNe la). 

The link of empirical relationships with the underlying explosion mechanism can be 
provided by approaches which, as the one presented here, involve density diagnostics: the 
late emission of SNe la "sees" through the ejecta at phases where spherical symmetry is still 
preserved. It is a good tracer of density profile and of central density of the ejecta, which 
can be plotted into WD mass (the central density/mass relationship is univocal for WDs). 
We thus expect that both the debate on the two alternative scenarios, Chandrasekhar or 
sub-Chandrasekhar, and the cosmic distance scale can be illuminated by these Co-powered 
nebular candles. 

Research on late emission of supernovae was started in collaboration with L.B.Lucy, 
whose advice on this and related subjects is highly appreciated. A large number of 
hydrodynamic models were explored thanks to the stimulating interest shown by Dave 
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Table 1: Distances from Nebular SNe la 

SN Galaxy Type (m-M)o E(B-V)» M^„^ Group / Cluster 

d (Mpc) 

SN 1989M NGC 4579 Sb 31.61± 0.2 0.0- 0.1 -19.11± 0.25 Virgo 

21 ± 2 

SN 1990N NGC 4639 Sb 31.81± 0.2 0.0- 0.1 -19.31± 0.25 Virgo 

23 ± 2 

SN 1992A^ NGC 1380 SO ^l-^toll 0.0-0.1 -19. 54^°;^^ Fornax 

24:11 

SN 1991T NGC 4527 Sb 30.57 ± 0.17 0.1-0.15 -19.53± 0.22 Virgo 

13 ±1 

SN 1991aa= Anl242-06 SB 33.12 ± 0.16 0.0-0.1 ...'^ vo=3300 " 

42 ± 3 

SN 1989B NGC 3627 SB 29.95 ± 0.4 0.3 -18.81 ± 0.45 Leo Group 

9.8 ± 1.5 

SN 1994D NGC 4526 SO 31.02 ± 0.25 0.0-0.1 -19.12 ± 0.3 Virgo 

16±2 



" Internal estimate of reddening 

^ This SN la might be below 1.38 M® 

= The coordinates of the galaxy are: a= 12^ 42" 34.7% (5=-06° 02' 38.3" (eq. 1950.0) 

unknown at maximum 
^ Vq: observed recession velocity of the parent galaxy in km s~^ 
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Figure Captions 

Figure 1 The electron density profile ol Chandrasekhar and sub-Chandrasekhar SNe la 
ejecta and its effect in producing different ratios ol Fe"'" forbidden emission. The [Fe II] 
A5262/A8617 ratio is displayed as a lunction ol ne, lor the maximum and minimum electron 
temperatures reached in those models (short dashed lines). For each model, the maximum 
electron density (ue) at each phase, here 270 (1) and 300 (2) days alter explosion) is given 
by the intersection ol the Ue profile (lines running Irom upper right to lower left) with the 
X-axis, and the projection ol this value on the curve ol the emissivity ratio gives an estimate 
ol the relative intensities ol the lines for that model. 

Figure 2 Comparison ol spectra ol SNe la (observed spectrum: solid lines; model spectra: 
dotted lines). The spectrum ol SN 1994D in NGC 4526 is plotted together with that ol the 
Chandrasekhar W7 model (Nomoto, Thielemann, & Yokoi 1984) at 300 days after explosion 
(top panel). A similar comparison is done (bottom panel) with the sub-Chandrasekhar 2D 
Model 6 Irom Livne & Arnett (1995), which yields an amount ol ^^Ni equal to that ol model 
W7 but in the explosion ol a 0.8 WD after accretion ol 0.17 ol He. Arrows indicate 
evidence ol a too low-density model Irom ratios ol forbidden line emission ol Fe"*". Stable Ni 
is deficient in sub-Chandrasekhar models, and it is, however, seen in "normal" SNe la. Its 
emission as [Ni II] A7379 is well accounted for by Chandrasekhar models. Among the class 
ol sub-Chandrasekhar models, model 6 is the one most closely resembling "normal" SNe la. 
It is representative ol He-detonations ol WDs ol total mass before explosion M final ~ 1^© 
which synthesize about 0.6 ol ^^Ni (equivalent models by other authors have similar 
features). Both less and more massive WDs exploding by He-detonations seem worse 
candidates for "normal" SNe la. 
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Fig. 1.- 
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SN 1994D in NGC 4526 
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